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O ■ ABSTRACT 

(N ; 

^ I Using the integration by parts method we derive a closed analytical expression 

/"V ' for the result of the integration of an arbitrary dimensionally regulated tadpole 

diagram composed of a massless propagator and two massive ones, each raised into 

' an arbitrary power, and including an arbitrary tensor numerator. We also briefly 

^ , discuss the implementation of the formula in the algebraic manipulation language 

^ ■ of FORM. 

o : 

Q . 1. Introduction 

(N : 

^ ' Feynman integrals (FI) witti complicated tensor numerators are usually dif- 

Q . ficult to work with. Even in the cases when they are known to be analytically 

^ i calculable in principle, in practical terms their evaluation often implies a tedious, 

P^l time-consuming and error-prone labour of reducing the problem to the calculation 

p • of a host of properly constructed scalar integrals. Moreover, the number of the 

O , latter integrals grows very fast as the numerator's structure gets more complicated. 

There exist only a few examples when the task is completely solved. A 
good example is provided with the so-called p-integrals, that is completely massless 
^ ■ Feynman integrals depending on only one external momentum. Here the explicit 

result is known for one-loop integrals^, while, say, a given three-loop tensor p- 
integral can be done by employing a rather cumbersome and time-costly method of 
harmonic projections^. 

In this talk we discuss another useful class of Feynman integrals — integrals 
without external momenta at all but comprising massive linesEi as well as massless 
ones. They will be referred to as m-integrals. Such integrals naturally appear in many 
problems where the mass m may be treated as a "heavy" one, much larger than all 
other mass scales involved. 

In the one-loop case m-integrals are rather trivial and we shall concentrate 
on two-loop m-integrals pictured in Fig. 1. 2-loop m-integrals with only one massive 
line (Fig. la ) may be reduced to 1-loop case after firstly integrating the 1-loop 
p-subintegral. 2-loop m-integrals with more than 1 massive lines (Fig. lb,c) are not 
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It is understood that all the massive propagators depend on one and the same mass m. 
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(a) (b) (c) 

Figure 1: Diffferent cases of 2-loop m-integrals: dashed lines are massless; solid lines have mass m. 

SO easy to do. We show how the use of the integration by parts method leads to a 
simple and general result for arbitrary (not necessarily scalar) two-loop m-integral 
with two massive and one massless line (see Fig. lb). In principle, our method 
allows also to reduce a tensor integral with 3 massive lines of Fig. Ic to a similar 
scalar integral. The latter can probably be done (at least in some particular cases) 
again through integration by parts*^. But to the best of our knowledge no explicit 
integration formula for this integral exists if the powers of all three propagators are 
arbitrary. 



2. Setting the problem 

We begin with a bit more complicated FI of the same topology as shown in 
Fig. lb but with a non-zero external momentum q. The corresponding analytical 
expression reads (in Euclidean space) 



A^h dPl2 V{p) 



(1) 



where £i and £2 are the loop momenta, p = {pi,P2,P3} are the propagator momenta 
and VQ is a tensor nominator. We shall deal with three possibilities of expressing 
the momenta p in terms of the loop momenta and the external momentum, viz., 

P2=^2, P3 = -(4+^2), (2) 
Pl^il, P2-^2+g, P3 = -(4+^2), (3) 

Pi^h, P2^h, P3 = q-{h+h)- (4) 
In a particular case of P{) = 1, g = the result for Eq. (|l]) is known"^ 

(Eq. ® with V{p) ^ 1, and q = 0) = {ni')D-»-P-y ^I^E^l^^^Mia, (3,j) 



with 



n ^ r(a + 7 - D/2)Ti(3 + 7 - D/2)Tia + (3 + ^ - D) 

M[a, P, 7 = — — — . 5 

r[a + fi + 2j - D) 

Our aim is to generalize this result on the case of arbitrary tensor polynomial 



Vip), still keeping q = 0. To simplify the formulas we shall set m = 1 below. 



2. Solution through the integration by parts method 

As is well known within dimensional regularization the value of the FI Eq. (|l|) 
does not depend on specifying the propagator momenta and all three choices 
are completely equivalent. The fact can be conveniently expressed through some 
differential identities. The essence of the integration by part method for dimension- 
ally regulated Feynman integrals consists of the use of such identities in order to 
simplify integrals to be computed^. 

Let us try to apply the method in our case. It is convenient to write eq. 
Eq. (|1|) in a condensed form as follows 

eq. Eq. (HD - / V{p)I{a,p,^) (6) 

and to introduce four differential operators acting on Eq. (^) 



Id . . 1 a 

29^ 



oTOi' * = 1,2,3 and 9^ - t^tttt- 



Explicitly, one has 



I(a B y) = — HRll — ll — M(a B y) 



(Pb)/(a, 7)) = (i^T^iP)) Ha, 7) - a,V{p)Iia + Su,P + + hi)pt , 

qf^ {V{p)I{a, P, 7)) = 1 ^P(p)^ I{a, /3, 7) - jV{p)I{a, /3, 7 + 1^, 

with ai = a,a2 = P, as = 7. The equivalence of the three momentum patterns (H-^) 
may now be expressed clS Si chain of identities 

p5'T(p)/(a,/3,7) = / pt^r{p)Iia,/3,j) = / /?, 7) - / 7). (8) 



These equalities allow us to evaluate immediately an integral of the form 
P.„»)/(o,^,7)l,.. - r(o)'rOTr(D/2) {(T)''^'-^''"} S(J^"<°'ft- + »>- 

with the Pochhammer symbol (a)„ = ^-feir^, -0 = 4-26 and V2n{p) being an arbitrary 

i [a) 

tensor in p = {pi,P2,p?.} of rank 2n. Indeed, in view of Eq. (H) we may freely replace 
Vinip) by V2n{q) = ^2n(pi = q,P2 = q,P3 = (?)■ Now the result Eq. @ comes from three 
simple observations: 

(i) The integral J I{a,p,j) is a scalar function of q and thus only the scalar 
component of the polynomial V2n{q) (that is proportional to g^") will survive after 
setting g = in the very end of the calculation. 

(ii) (^)"(g)"=n!(2-e)„. 



(iii) {^y^ = ilUl - 1 + e)« = (-)"(7)n(2 - 6 - 7 - n),,j^^. 

Thus we are left with the task of finding a representation of the initial integral 
Eq. (|l|) as a linear combinations of integrals of the form displayed in Eq. (^). The 
problem is solved by the use of the following identity 

["/2] 

rn{q)f{q') = E Vir {°^^"(«~)} f^-^^^Hi') (10) 

a=0 

where f{x) is an arbitrary smooth function of x and /'■"■^x) is defined in such a way 
that 

(11) 

It should be clear now that once the initial integral Eq. (|l]) has been expressed as 
a linear combinations of integrals of the form Eq. (P) it may be done without any 
problem. Indeed, without essential loss of generality we may assume that the poly- 
nomial V{p) = "Pni.nsCp) does not depend on and meets the following homogeneity 
equation: 

Now a direct application of Eq. , Eq. (H) and Eq. ( p!0| ) gives 

y 'Pni,n2{Pl,P2)I{a,(3,-f) = 
rfJl/O \ [("i+»2)/2] [n2/2] [ni/2] 

C_')("l+"2 + 0'3) 

,(;,+.2+.3) , {°P2^°p^^n..n.bl,P2)} U=p2=P3} Ip3=0 

which is the formula we wanted. 

The algebraic structure of Eq. (|T2|) is very similar to that of the corresponding 
formula for 1-loop j3-integrals in^. This observation has allowed us to perform a 
simple algebraic programming of Eq. (|12|) in FORM^ by closely following the routine 
ONE.PRC from the package MINCERA 
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